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PACS 74 . 20 . Mn - Nonconventional mechanisms 

PACS 74 . 20 . -z - Theories and models of superconducting state 

PACS 71.27.+a - Strongly correlated electron systems; heavy fermions 

Abstract -We perform a renormalization group (RG) study up to two-loop order of an effective 
low-energy two-band model to describe some of the recently discovered iron-based superconduc- 
tors. Our starting point is the itinerant electronic model proposed by Chubukov et al. [Phys. Rev. 
B 78, 134512 (2008)], which displays two small, almost nested Fermi pockets with one hole pocket 
centered at (0,0) and one electron pocket centered at Q = (n,n) in the folded Brillouin zone. 
We then proceed to implement a complete two-loop RG calculation for this model of four-point 
vertex corrections, quasiparticle weight and several order-parameter susceptibilities in order to 
evaluate the robustness of one-loop RG results available in the literature with respect to including 
self-energy effects and higher-order quantum fluctuations. 
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Introduction. - The recent experimental observation 
2] that the iron-based pnictides (such as LaFeAsO and 
SrFe2As2, to name but a few) exhibit unconventional su- 
perconductivity at critical temperatures up to T c = 55K 
sparked big efforts in the community of strongly corre- 
lated systems to understand these materials. This is be- 
cause they are the first high-T c superconductors ever dis- 
covered which are outside the well-known copper-oxide 
family. Since electron-phonon coupling in these iron com- 
pounds appears to be very small [3], a purely electronic 
mechanism for superconductivity emerges as a strong pos- 
sibility. Indeed, in many aspects, the pnictides resem- 
ble the physics of the cuprate superconductors: these two 
classes of compounds are effectively two-dimensional ma- 
terials which exhibit antiferromagnetism U at zero dop- 
ing and both become a superconductor upon doping. This 
observation could imply that the discovery of the underly- 
ing mechanism of superconductivity in the iron pnictides 
might also give important insight into solving the long- 
standing cuprate high-T c superconductivity problem. 

By contrast, there are some crucial differences between 
the pnictides and the cuprates, which makes the former 
materials also very interesting from a fundamental point 
of view. Unlike the cuprates which always display a local- 
ized Mott insulating phase at low doping, the pnictides are 
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instead semi- metals for this doping regime [3]. Besides, 
there is growing consensus in the community that some 
of these iron-based materials remain itinerant for all dop- 
ing levels [5]. This suggests that the pnictides are in fact 
less correlated than the cuprates and, therefore, weak cou- 
pling theories could be a valid starting point to describe 
at least qualitatively the properties of these former mate- 
rials [H[2]. Moreover, from angle-resolved photoemission 
spectroscopy (ARPES) experiments [HIH] and theoretical 
band-structure considerations [TOJ[TT], it is now generally 
accepted that the low-energy electronic structure of these 
materials are in some respects similar to one another, with 
a Fermi surface consisting essentially of two small electron 
pockets centered around M = (tt, 7r)-point and two (or 
three) small hole pockets centered around r = (0, 0)-point 
in the so-called folded Brillouin zone. 

Motivated by these experimental results, many re- 
searchers have put forward various types of multi-band 
electronic models aiming to explain some of the key prop- 
erties displayed by these materials, most notably, their 
phase diagram and also the symmetry of the corresponding 
superconducting gap. These studies included two-band 
[12], four-band [IMS] and five-band models [17]. These 
models typically contain several competing ordering ten- 
dencies at low-energies, which must be in principle treated 
on equal footing. In this respect, a promising theoretical 
framework which is tailored for describing such systems 
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Fig. 1: (Color online) The Fermi pockets of the two-band model 
analyzed in this work. We also show schematically the inter- 
action processes included in this model: Ui and V% stand for, 
respectively, interband couplings with (0,0) and Q = (tt, tt) mo- 
mentum transfer, U3 corresponds to interband pair scattering, 
and U4 refers to all intraband interactions. 

turns out to be the renormalization group (RG) approach 
in view of its unbiased nature |18j . 

Soon after the discovery of the pnictides, Chubukov et 
al. proposed a simple but ingenious effective two-band 
model [12] , which displays two small, almost nested Fermi 
pockets in the folded Brillouin zone to describe such ma- 
terials. By applying a one-loop RG approach - which is 
essentially equivalent to summing the so-called parquet 
diagrams up to infinite order - they concluded that at 
zero doping the model displays a antiferromagnetic spin- 
density wave (SDW) phase in agreement with experiments 
and, upon doping, a superconducting (SC) phase emerges 
with a extended s^-wave gap symmetry, as first suggested 
in Rcf. . In another important work [17j done by Wang 
et al., a more complex five-band model for the iron pnic- 
tides was analyzed using a functional generalization of the 
one-loop RG approach, and their results also point to an 
antiferromagnetic phase for the undoped system and a ro- 
bust s ± -wave pairing state which arises at larger doping. 
In order to further demonstrate that these results did not 
depend critically on some details of the band-structure of 
the above models and the approximations used, the au- 
thors in Rcf. [T5] set out to discuss a four- band model 
within the one-loop functional RG treatment and their 
data provided yet another confirmation of the previous 
RG results regarding SDW and s ± -wave symmetry super- 
conductivity as leading instabilities obtained within the 
two-band and five-band models. 

Another aspect of these iron-based superconductors 
concerns their gap structure i.e. the k-dependent varia- 
tion of the SC order parameter within a given symmetry 
class - and the related question of existence (or absence) 
of nodes on them. Two works recently addressed this is- 
sue within a one-loop RG approach by comparing the re- 




Fig. 2: Sunset diagram for the self-energy at two-loop order 
which yields a non-analytic contribution for energies larger 
than the Fermi energy in the present two-band model. 

suits obtained within a four-band and a five-band model 
[THHO]. As a result, Thomale et al. concluded that in 
the four-band case there should appear nodal gaps on the 
electron pockets, while the hole pockets are always node- 
less. By contrast, if an additional hole pocket is included 
in the model (hence, resulting in a five-band model), the 
net effect of this would be turning the nodal gap on the 
electron pockets nodclcss. This could explain why some 
compounds seem to exhibit experimentally nodal gaps and 
others not. This conclusion is shared by the authors in 
Rcf. [20] who also gave a thorough analysis of the two- 
band model case. In this latter work, it was shown that the 
two-band model naturally describes nodeless gaps on both 
electron and hole pockets which agrees qualitatively with 
the five-band scenario. This could suggest that both two- 
band and five-band models might be in the same universal- 
ity class and therefore have the same low-energy physics. 
Consequently, this would imply that the two-band model 
might be a good, minimal low-energy effective model in 
order to describe some iron-based superconductors which 
display nodclcss gaps on all pockets. 

For this reason, we revisit in this work the two-band 
model discussed previously by Chubukov et al. within a 
one- loop RG approach. As a result, we report here a com- 
plete two-loop RG calculation for this model of all vertex 
corrections, self-energy and several order-parameter sus- 
ceptibilities with the important goal of evaluating the ro- 
bustness of one-loop RG results available in the literature 
with respect to including self-energy effects and higher- 
order quantum fluctuations. 

Model. - Our starting point is the symmetric phase 
of the two-band model proposed in Rcf. [T2] that displays 
two small, almost nested Fermi pockets - with one hole 
pocket centered at T-point and one electron pocket cen- 
tered at M-point in the folded Brillouin zone (see Fig. 1) - 
and which includes also both interorbital and intraorbital 
interactions. If we use a coherent-state functional integral 
representation of the resulting Hamiltonian after a suit- 
able canonical transformation (for a thorough discussion 
on this point, see Ref. [E]), the model at T = and con- 
stant chemical potential /i = Ep becomes described by 
the action S = Sq + Si n t , where 

So = ^2[{-iuj + e c p )c pa c pa + (-iu + e P )J pa f pa ], (1) 

p,cr 

and 
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where P4 = Pi + P2 — P3 and the volume V has been 
set equal to unity. We assume here, for simplicity, that 



p 2 /2m and e 



/ 

p+Q 



-e£. Besides, 



prr 



and 



c pcr are, respectively, the creation and annihilation Grass- 
mann fields for fcrmions with spin projection a that are 
near k = (0, 0), f pa and / p(T are the creation and annihi- 
lation Grassmann fields for fcrmions with spin projection 
a that are close to Q = (7r,7r), and U> (for i = 1..4) are 
the microscopic (bare) coupling constants of the model 
which are displayed schematically in Fig. 1. The above 
action defines our bare quantum field theory which must 
be regularized in the ultraviolet by restricting the energies 
to |w| < A (i.e., W = 2A represents the bandwidth of 
the model). 

Since both hole and electron pockets of the iron pnic- 
tides are small compared to the bandwidth W as measured 
from ARPES experiments, we shall concentrate through- 
out this work only on the physical regime of the two-band 
model in which the energies are actually larger than the 
Fermi energy (\uj\ > Ep). In this way, if we calculate 
the self-energy of the model up to two-loop order for this 
regime, we obtain that its non-analytic contribution (Fig. 
2) is given approximately by 



E(iw, p = kp) 



(0)l2 



A^(oWln( ^ 



(3) 



with kp being the Fermi vector and N(0) = m/2ir is the 
constant density of states of a two-dimensional Fermi gas. 
The presence of this non-analyticity in the self-energy is 
a generic feature of this model and several logarithmic 
divergences also appear if one calculates four-point ver- 
tex corrections and order parameter susceptibilities within 
perturbation theory up to two-loop order. 

Method. - The field-theoretical RG approach up to 
two loops that we shall use now is standard [21] and 
details also appear elsewhere in the context of another 
fermionic model |22j . In order to circumvent the problem 
of logarithmic singularities and non-analyticities emerg- 
ing in the low-energy limit of the present model within 




Fig. 3: Schematic representation of the Feynman diagrams of 
vertex corrections included in the present two-loop RG calcula- 
tion. The single-particle propagators are represented by either 
solid or dashed lines according to their association with the cor- 
responding Fermi pockets. The diagrams with crossed squares 
represent the counterterms. 



perturbation theory (see Fig. 3), the field-theoretical RG 
strategy consists of rewriting the unobserved bare quan- 
tities of the microscopic model defined by Eqs. (JlJ and 
([2jl in terms of the experimentally observed renormal- 
ized parameters plus appropriate counterterms. These 
counterterms have the main effect of regularizing classes 
of diagrams to a given order in the model at a float- 
ing new energy scale A, such that the rcnormalizcd per- 
turbation theory becomes well-defined in the low-energy 
limit (i.e. A — > 0). An important point we wish to 
stress here is that both coupling constants and fermionic 
fields of the model must be renormalizcd at two-loop 
RG level. Therefore, we must define: c pcr = Z 1 / 2 Cp Cr , 

f P a = ^ 1/2 /pl, and U^ 0) = N-^tyZ- 2 ^ + A Ul ), where 
Z = (1 — dT,(iuj,p = kF)/<9(iw)| w= o) -1 is the quasiparti- 
cle weight, Cp CT and fp a stand for the renormalized fields, 
Ui (for i = 1..4) represent the corresponding dimensionless 
renomalized coupling constants of the model, and Aui re- 
fer to the counterterms which must be calculated order by 
order in perturbation theory. Since this program is suc- 
cessfully accomplished, the field theory model analyzed 
here is indeed renormalizable. 
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Results. - We can now adjust the counterterms Aitj 
such that all divergences are exactly canceled in our series 
expansion up to two-loop order. But the price we pay for 
this is the appearance of a new scale A with all physical 
quantities now depending on this scale. By contrast, the 
original model has no information about this quantity, i.e. 
the bare parameters do not depend on A. This leads us 
to the renormalization group conditions for the bare cou- 
plings of the model, i.e. A(dU^ / dA) = 0. As a result, 
the RG flow equations for the renormalized dimensionless 
couplings at two-loop order become 



1*1 



+ {u\ + uf) + ( u 

— 2u±(u\ + 1*2 — - 



L - U 2 )u\ 

u\u 2 + -uf), 



"2 

"3 = 
«4 



2(ui - u 2 )(u 2 + u 2 ua) - 2u 2 u\, 
= + 2u 3 (2ui - u 2 - it 4 ) - 1*3, 

(«3 + «2) 



(4) 

(5) 
(6) 



114U3 



2u 2 us 



2u\ 



+ 2uii* 3 + 3i*2i* 3 — 2u\ — 2i*ii* 2 + 2u l u 2 , (7) 

where the derivatives are taken with respect to £ = 
(1/2) ln(Ao/A). The initial conditions for this system 
of differential equations are naturally given by the mi- 
croscopic interactions, i.e. i*;(£ = 0) = N(0) (for 
i = 1..4). A schematic representation of the Feynman 
diagrams corresponding to the vertex corrections up to 
two-loop order are displayed in Fig. 3. 

First, we focus our attention on the numerical solution 
of the flow equations for the renormalized couplings as a 
function of A. This was performed by means of fourth- 
order Runge-Kutta method. We analyze simultaneously 
both one-loop and two-loop RG flows for a direct compar- 
ison of the two approaches. In both cases, even though all 
couplings indeed diverge at low energies (£ — > 00), their ra- 
tio approach infrared (IR) stable fixed points in this limit 
(Fig. 4). Most importantly, as can be clearly see from 
this figure, the inclusion of two-loop order fluctuations in 
the RG scheme has the effect of changing the fixed point 
structure of the field theory model. 

In the one-loop RG approach, as first obtained in 
Ref. |12j . there are two IR stable fixed points given by 
(1*2/1*1) — > 0, 1*3 = ±\/5i*i, 1*4 = —iti, which control the 
low-energy dynamics of the model. Therefore, from the 
RG analysis of the corresponding susceptibilities up to one 
loop, the authors in Ref. [T2] concluded that the leading 
coupling in this limit turns out to be the 1*3 pair scatter- 
ing which appears to be the main responsible for inducing 
a SDW instability at zero doping and which gives rise to 
an extended s ± -wave superconducting instability at larger 
doping in the model. 

On the other hand, in the complete two-loop order RG 
approach implemented here in this work, there is only 
a single IR stable fixed point which is given instead by 
(1*2/1*1) — > 0, (1*3/1*1) — > 0, 1*4 = —i*i. Hence, higher- 
order fluctuations shift the dominance of the coupling con- 
stants. Although the 1*3 coupling still diverges in the low- 




Fig. 4: (Color online) RG flow at one-loop order (upper panel) 
and two- loop order (lower panel). (Note that the axes in the 
above plots are different.) 



energy limit, it does so at a slower rate compared to the 
one-loop RG approach. As a consequence, the interaction 
processes given by the interband and intraband forward 
scatterings (i.e. u\ and 1*4, respectively) eventually over- 
come the 1*3 pair interaction in this limit. It is interesting 
to note that the new IR fixed point found here at two- 
loop RG level still satisfies the SO(6) symmetry condition 
|23j . Moreover, the quasiparticle weight Z renormalizes 
very weakly in this regime and always remains closer to 
unity, thereby indicating Fermi liquid behavior. 

In order to verify if the interband and intraband forward 
interactions are able to induce density-wave and pairing 
instabilities in the model, it is important to calculate the 
corresponding susceptibilities by introducing an infinites- 
imal external field in the appropriate channel and evalu- 
ating its linear response. Therefore, we must add to the 
action that describes the present model the following term 

Sext = ^""^ Tg C c Ck,g-C_k,-o- + ^""^ TsC. f /k,q/-k,-cr 



c k,a/k+Q,/j + h.C. 



(8) 



where and T^^ P are the bare response vertices 

for the superconducting and density-wave orders, respec- 
tively. This added term will generate new Feynman di- 
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Fig. 5: Feynman diagrams for the response vertices for both 



superconducting [T^q c ^\] and density-wave [7^/*"] orders. 
The diagrams with crosses represent the corresponding coun- 
terterms. 



agrams - the three-legged vertices displayed in Fig. 5 - 
which will also generate new logarithmic singularities in 
the low-energy limit of our field theory model. Therefore, 
we must regularize these divergences (see, e.g., Ref. |24j 
in the context of another fermionic model) by defining 
the rcnormalizcd response vertices and the corrcspond- 



-(0) 



= Z~ 



l (Tsc,c(f) 



ing counterterms as follows: Tg C c ^ 

AT S o,c(/)) and T { X\~ Z~\T^w + A7^). Again, by 
invoking the RG condition for the bare quantities of the 
model A(dT- 0) /dA) = (for i = SC and DW), we obtain 
the RG flow equations for the response vertices 



T s 



C,c(/) 



u 4 



y ) Tsc,c(f) " u 3T~scj(c), 



1 DW 



"1 



DW 



-U 3 Sa-piTjDw]*, 
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with a,f3 =t, I, and the derivatives here are also taken 
with respect to £. By symmetrizing these response ver- 
tices, we obtain the following order parameters 



1 SC 



Tsc,c — Tscj 
= Tsc,c + Tsc,f, 

T - -rtt 

I CDW(SDW) ~ 

7~CDW{SDW)± 



T'YX _i_ 
1 DW ^ 'DW' 



TcDW(SDW) ±^CW(SDW)' 



where Tg S c ' and T^q represent, respectively, extended 
s-wave and conventional s-wave superconducting orders, 
Tcdw and Tsdw stand for charge density wave and spin 
density wave correlations (the subscripts + and — refer 
to, respectively, the real and imaginary parts up to a con- 
stant of the corresponding order parameter). The RG flow 
equations at two loops then become 




Xsdw+ 



two-loop RG 
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Fig. 6: (Color online) RG flows of several susceptibilities of the 
model for both one-loop (left panel) and two-loop (right panel) 



RG approaches. The initial conditions used were uf^ — u^ 1 ' — 
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(12) 

2u 2 T u 3 - yj Tcdw±, (13) 

(14) 



T (0) for all 



«3 + «4 + y )T [ s s I 



Ul 



Ml ± U 3 - — 75DW± 



0) 



with initial conditions given by 7i(£ 
order parameters defined above. 

Once we computed the response vertices associated with 
an instability towards a given ordered phase, we can then 
proceed to calculate their corresponding susceptibilities. 
They are given by 

Xi(0 = f dtTi(0V(0, (is) 



where % = SC - s ± , SC - s, CDW± and SDW±. In Fig. 
6, we plot the results both for one-loop (which agrees with 
Ref. [2D]) and the two-loop RG approach implemented 
here. We observe that the instabilities of the model within 
our RG scheme are surprisingly not changed qualitatively 
compared to one-loop results, even though the divergence 
at two loops takes place more slowly from a quantitative 
point of view. 

Conclusion. - We have carried out a complete two- 
loop RG study of an effective low-energy two-band model 
to describe some recently discovered iron-based supercon- 
ductors. We have shown that the inclusion of two- loop 
quantum fluctuations has the main effect of changing the 
fixed point structure of the model. Even though the main 
instabilities turn out to be qualitatively the same as in 
one-loop calculations, the present work suggests a different 
microscopic mechanism with the dominant interactions at 



p-5 



Vanuildo S. de Carvalhol Hermann Frcirel 3 



low energies being the intraband and interband forward 
interactions. It would be very interesting to implement 
such a two-loop RG scheme also at finite temperatures in 
order to calculate other physical quantities such as the uni- 
form charge and spin susceptibilities of the model and to 
further compare the theoretical predictions obtained from 
the present RG approach with recent experimental studies 
[251126] performed on these important materials. 
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